The exact analytical solution of the linear structure growth rate in ACDM cosmology 

and its cosmological applications 



Pengjie Zhang^Q 

^ Key laboratory for research in galaxies and cosmology, Shanghai Astronomical Observatory, 
Chinese Academy of Science, 80 Nandan Road, Shanghai, China, 200030 

We derive the exact analytical solution of the linear structure growth rate in ACDM cosmology 
with flat or curved geometry, under the Newtonian gauge. Unlike the well known solution under the 
Newtonian limit (Heath 1977, our solution takes all general relativistic corrections into account 
and is hence valid at both the sub- and super-horizon scales. With this exact solution, we evaluate 
cosmological impacts induced by these relativistic corrections. (1) General relativistic corrections 
alter the density growth from z — 100 to 2 = by 10% at fc = 0.01/i/Mpc and the impact becomes 
stronger toward larger scales. We caution the readers that the overdensity is not gauge invariant 
and the above statement is restrained to the Newtonian gauge. (2) Relativistic corrections introduce 
a scale dependence in the density fluctuation. It mimics a primordial non-Gaussianity of the 
local type with /nl''' ~ 1. This systematical error may become non-negligible for future all sky deep 
galaxy surveys. (3) Cosmological simulations with box size greater than IGpc are also affected by 
these relativistic corrections. We provide a post-processing recipe to correct for these effects. (4) 
These relativistic corrections affect the redshift distortion. However, at redshifts and scales relevant 
to redshift distortion measurements, such effect is negligible. 
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INTRODUCTION 



tor k satisfies 



A well known and widely used solution of the linear 
density growth in ACDM cosmology, first derived by 
Heath 1977 [ij, is 



, oc Dm,N oc H 

Jo 



da 



(1) 



Here, 5m is the matter overdensity. Dm,N is the linear 
density growth rate in the Newtonian limit. Throughout 
the paper, the subscript "N" denotes the corresponding 
property in the Newtonian limit. H — H(a) is the Hub- 
ble parameter at redshift z = 1/a — 1 and a is the scale 
factor. This solution is derived under the limit of negli- 
gible radiation, negligible baryon pressure and identical 
initial conditions for fiuctuations in baryons and dark 
matter. These conditions are adopted throughout the 
paper. 

This solution is valid for arbitrary cosmological con- 
stant A. It is also valid for arbitrary curvature K, at 
scales much smaller than the curvature radius tk — 
= 1/ HQ^J\^lK\■ Here is the dimensionless 
curvature density and i/o = H{a = 1) is the present 
day Hubble constant. Throughout the paper, we set the 
speed of light c = 1, so 1/Hq — 3/i^^Gpc. 

However, this solution is derived by neglecting rela- 
tivistic corrections to the Poisson equation and to the 
continuity equation, so it is valid only in the Newtonian 
limit and hence at sub-horizon scale, where the wavevec- 



3 X lO^/i-iMpc 



(2) 



'Electronic address: l pjzhang@shao.ac.cn| 



Here, fig and JIa are the dimensionless matter and A 
today, respectively. 

Modern surveys are pushing the observational bound- 
ary to larger scale and higher redshift. This enables 
several important cosmological applications. (1) Galaxy 
clustering at scales k ^ 0.01/i/Mpc is a sensitive probe 
of primordial non-Gaussianity (PMG). It relies on the 
asymptotic behavior of galaxy clustering at A: — to 
probe PMG (e.g. 0). (2) The large scale structure (LSS) 
growth history, often inferred by combining low redshift 
and high redshift measurements, contains key informa- 
tion to test general relativity and probe dark energy. 
Such measurement is usually done at k ^ 0.01/i/Mpc 
scale. However, since the horizon scale decreases with 
increasing redshift, even a mode deep in the sub-horizon 
regime today can be close to the horizon scale at earlier 
epoch. 

For these cosmological applications, whether relativis- 
tic corrections can be safely neglected and whether Eq. 
[1] is sufficiently accurate are becoming issues of practi- 
cal importance and have evoked discussions (e.g. [3|[). 
Existing packages such as CMBFAST and GAME [1 
numerically solve the fully general relativistic linear per- 
turbation equations. So they can be reliably resorted to 
calculate the structure growth at all linear scales. 

However, it is still useful to understand the role of these 
relativistic corrections through analytical approach. This 
motivates us to derive the exact solution of the linear den- 
sity evolution with all relativistic corrections included. 
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valid at both the sub-horizon and super-horizon, for ar- 
bitrary A and K . 

For LSS study, the most widely used gauge is the New- 
tonian gauge (or the conformal Newtonian gauge) . Hence 
we will perform the derivation of the exact solution in 
the Newtonian gauge, and then discuss its implications 
on the LSS-based cosmology. Some of the perturbation 
variables in the Newtonian gauge, such as the matter 
overdensity (5™, are not gauge invariant. However, such 
gauge freedom is not a problem, as long as we prop- 
erly connect the perturbation variables in the Newtonian 
gauge to observables (refer to [6| for a comprehensive 
treatment). Furthermore, once the solution in the New- 
tonian gauge is derived, it can be transformed to other 
auges by performing a coordinate transformation (e.g. 
I). For these reasons, we will stick to the Newtonian 
gauge throughout the paper. 

After we derived the solution and were writing this 
manuscript, we found an earlier work by Jai-Chan Hwang 
0, who has derived equivalent solutions in not only the 
Newtonian gauge, but also in other five gauges. Nev- 
ertheless, we feel the need to publish our result, for a 
number of reasons. First, the solutions are derived in 
independent ways. [1] first derived the solution in the 
comoving gauge and then converted to the Newtonian 
gauge (denoted as the zero shear gauge in [§]). For us, 
we derive the solution directly in the Newtonian gauge. 
Our result thus serves as an independent check and the 
agreement between the two verifies the validity of both 
solutions. Second, given recent discussions on the role 
of these relativistic corrections (e.g. [1]), it is useful to 
address the existence of these exact solutions and clar- 
ify some confusions arisen in the literature. Third, the 
exact solution derived has direct applications to modern 
cosmological topics, some of them, e.g. the primordial 
non-Gaussianity, will be elaborated later in the paper. 

The exact analytical solution that we found is valid 
for any A and curvature K . But numerical results and 
related discussions are based on the standard flat ACDM 
with fio = 0.268, IIa = 1 — ilo and Q,k = throughout 
the paper, unless otherwise specified. 



will be used later, 



K = inGpma^ + HH a 
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(5) 



Throughout the paper, = d/da. 

For reviews of linear perturbation in ACDM cosmol- 
ogy, please refer to, e.g, 7, 8, 10]. Here we just summa- 
rize relevant results. For the convenience of LSS study, 
we adopt the Newtonian gauge defined through 

ds^ = -(1 + 2(t))dt^ + a^(l - 2(t))-/^jdx'dx^ . (6) 

Here, we have adopted the result that the two Newtonian 
potentials are identical in ACDM (e.g. Q). We choose 
the following three independent perturbation equations. 



2 tt2, 



aH{(j) a + 4))^ 'i-KGpmO^W , 



(7) 



(8) 



(9) 



Here, = 7'-' V;iV;j and the covariant derivative V;i is 
defined in the 3D constant curvature space, v = —VW 
is the peculiar velocity and W is the velocity potential. 
The three equations above compose a complete set for the 
evolution of the three perturbation variables, the New- 
tonian potential cj), the matter overdensity 5m and the 
velocity potential W . 

The (f) a + (f) term in Eq. [7] is a relativistic correction 
to the Poisson equation. The cj) term in Eq. [9] reflect 
the ambiguity to distinguish inhomogeneities in space- 
time metric to that in the matter-energy fluid. These 
relativistic corrections are usually neglected in LSS study. 
It is the goal of this paper to quantify their impacts. 



Linear evolution of the gravitational potential and 
peculiar velocity 

Interestingly, including all these relativistic terms al- 
lows us to derive the exact solution. From Eq. [71 we 
obtain , 



THE EXACT SOLUTION 

We work with the FRW background metric 

ds^ = —g^^dx^dx'^ = —dt^ + a^^ijdx^dx^ 
dr^ 



-dt^ 



'de^ 



1 - Kr^ 

The Hubble expansion rate is given by 

rr2 SttG , ^ K 
H = (Pm+PA) ^ • 



(3) 



(4) 



Here p,„ is the mean cosmological matter density. Taking 
the derivative of both sides, we obtain a relation which 



(V2 + 3K)(j) - 3a^H^{(j) a + 0) 



(10) 



AirGpma'^ 

Plug the above equation and Eq. [5] into Eq. [HI we obtain 



(V2 + 3K)(I) ~ 3a^H^{(j) a + 
AirGprna^ 



a + (j)) 



inGpma^ 



f 3(/. 
(11) 

Since pm oc a ^, in the above equation, the two terms 
cx cancel exactly. We then have 



3K(j)-3a^H^{(l) a + tjj) 
AirGprna'^ 
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(12) 
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The solution is 



a + (t>)~ K(t> 



— constant . 



(13) 



Immediately we see the linear evolution in is scale inde- 
pendent. This tells us that the solution of (j) obtained in 
the sub-horizon limit is also valid at super-horizon scale. 

Multiple both sides by inGpma^ / H^a^ and use Eq. O 
we obtain 



a(j) 
ll 



1 



(14) 



Here, we have used the fact that ATTGpmO?' is a constant. 

The above equation has two independent solutions. 
The solution (p oc H/a corresponds to the decay mode 
in density. So it is of little practical interest and will not 
be discussed hereafter. The other solution corresponds 
to the growth mode in the density evolution, 



H r da 

' cx cx — 



a Jo H^a? ■ 



(15) 



Now the scale independence in is explicitly shown. 

The linear evolution in velocity (divergence) can be 
obtained combining Eq. [5] & [T5l 



Ha 1/H^a^ 
IT ^ ^^da/H^a^ 



Here, 



H a 



1/H^a^ 



dill ( 



H rda/ma^ d\n< 



da 

■ 

(16) 
(17) 



For the standard ACDM, this quantity approaches unity 
at high rcdshift and decreases to ~ 0.5 at z = 0. 

Wc notice that both the evolution in the gravitational 
potential and in the velocity are identical to the ones 
derived in the sub- horizon (Newtonian) limit, namely, 



D,, = D. 



<I>,N 



D,, = D 



v,N 



(18) 
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FIG. 1: The relativistic correction C'{k,a). Relativistic cor- 
rections alter the matter power spectrum by a factor (1-1- C)^. 
So a positive C means the amplification of the matter cluster- 
ing. Since C oc we only plot C{k, a) at k — 0.001/i/Mpc. 
At 2; ^ 1, C{k, a) oc since dln(j>/dlna — >■ 0. 



Only when the scale is much smaller than the curvature 
radius and the horizon scale, we recover Eq. [1] 

CMB observations, along with other probes, show the 
universe to be nearly flat 



\^k\ < 1) [Hi- Hence aU 
modes accessible to observations are much smaller than 
the curvature radius. We then neglect the curvature term 
and replace with the usual Laplace operator in the 3D 
Euclid space. We further proceed to the Fourier space 
(V^ — >■ — fc^). We denote the corresponding Fourier com- 
ponents with a superscript . We then have 



k^4>x [1 + C(fc,a)] 



(20) 



Indeed, this is what expected from the scale indepen- 
dence in Dri, and D„. 



The impact of relativistic corrections is then completely 
captured by the relativistic correction term 



Linear evolution of the matter overdensity 

However, the situation for the overdensity evolution is 
different. From Eq. [71 



dm = -rrrrrr X 



(19) 



3K - 3a^H 



2tj2 



H a 



l/H^a' 



H rda/H^a^ 



C{k,a) 



3a^H^ Ha 



l/H^a^ 



fc2 y H da/H^a^ j 

a^{H/Hof (Ha l/H^a^ 



(21) 



3{k X 103/i-iMpc)2 \^ H J^da/H^a^ J ' 

For the standard ACDM cosmology, C > 0. The numer- 
ical result is plotted in Fig. [TJ 

The matter power spectrum Pm is related to the power 
spectrum of the gravitational potential through Eq. 
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FIG. 2: The linear density evolution in ACDM under the 
Newtonian gauge, for various k. The solid line corresponds 
to the one in the Newtonian limit. Larger deviation from 
the Newtonian result happen at higher redshifts, where the 
relativistic corrections are larger due to smaller horizon scale. 



k^P4>x[l + C{k,a)f (22) 



= Prn,N{k,z) X [l + C{k,a)] 



Since C > 0, relativistic corrections enhance the matter 
clustering. 

The linear density evolution is derived combining Eq. 



H 



da 



X [1 + Cik,a)] 



(23) 



Relativistic corrections induce a characteristic scale de- 
pendence fc~^ in the otherwise scale independent density 
growth in ACDM cosmology. Only when k ^ aH, does 
this scale dependence vanishes {C{k,a) <^ 1). The lin- 
ear density growth factor Dm is often defined as Dm = 
Sm{o-)/Sm{a = 1) and normalized at z = 0, in LSS study. 
We then have 



Dm (a, k) 



H{a) 



j;da/Hh 



Ho 

DmM{a) X 



1 + C{k,a) 
1 + C{k,a= 1) 



1 + C{k,a) 



1 + C\k,a= 1) 



(24) 



FIG. 3: Same at Fig. [2] but only show the evolution at z < 10 
where LSS surveys may have access. 



We plot D{a, k)/a for various k in Fig. [2] & [3] . As 
expected, the relativistic corrections become bigger at 
higher redshifts and can completely invalidate the New- 
tonian result (Eq. [1]), due to shrinking horizon. This an- 
alytical result verifies the finding of [3] , who used numer- 
ical calculation to quantify the relativistic corrections. 
One thing to notice is that, [3] normalizes the overdensity 
at z = 1100 and hence defines D = 5m{z) / 5m{z = 1100), 
so deviations from the Newtonian limit appear at lower 
redshift instead. 

Surprisingly, the large deviation shown in Eq. [T] and 
Fig. [5] does not necessarily invalidate the applicability of 
the Newtonian approximation (Eq. [IJ in LSS study. This 
is essentially a normalization issue. To better demon- 
strate this point, we highlight the linear density growth 
at z < 10 in Fig. [31 Most LSS surveys are limited to 
z < 4. If we solely compare LSS at redshifts accessible 
to those surveys to measure the structure growth rate, 
Eq. [l]is essentially exact for k > 0.01/i/Mpc. Even for 
the mode k — 0.006ft./Mpc, the accuracy is better than 
1%, negligible comparing to the cosmic variance. Hence 
in general, Ec^. [T]is sufficiently accurate to describe the 
linear growth at redshifts accessible to LSS surveys at 
z ^ 4. However, later in the paper we will discuss an im- 
portant exception, where by the multiple tracer technique 
proposed by [l^ . the cosmic variance can be eliminated. 

21 cm surveys can probe LSS to the reionization epoch 
at z 10 and even higher [l3| . A comparison between 
these epochs and z ~ based on the Newtonian approx- 
imation (Eq. [1]) would lead to ^ 1% error in the derived 
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FIG. 4: General relativistic corrections to the Poisson equa- 
tion and the continuity equation can mimic a primordial 
non-Gaussianity of local type and induce a bias oc 
bg/(bg ~ 1). The plot shows the case for bg = 2 and bg = 0.5 
respectively. 



Primordial non-Gaussianity 



The galaxy overdensity dg is an observable and hence 



gauge independent, 
expressed as 



In the Newtonian gauge, it can be 



bgS„ 



(25) 



Here, / is a linear function of (f> and v. We have shown 
that the relativistic corrections discussed earlier do not 
bias the Newtonian calculation on (/) and v, so the only 
term affected by the relativistic corrections is the usual 
bgSm term, where bg is the galaxy bias. 

Primordial non-Gaussianity (PMG) of the local 
type induces scale dependence to the otherwise scale- 
independent bias at large scale [3], 



Abg^2{bg^l)f^°r'''5c ™ " 



" 2ag{a)k'^ 



(26) 



Here, /^l^' is the PMG parameter. 5c is the threshold 
for halo collapse and g oc D(k — ^ oo)/a, normalized such 
that g{a — > 0) — > 1. If we miss the relativistic corrections 
in the matter evolution, it will bias the interpretation of 
galaxy bias by 



6bg — bgC{k, a) . 
This in turn causes a systematical error in /^l^', 



(27) 



density growth rate at fc = 0.01/i/Mpc, still small, but 
may no longer be negligible. 

On the other hand, if we want to compare LSS mea- 
surements at low redshift to CMB at z ~ 1100, Eq. 
[T] is no longer applicable, shown in Fig. [5J Even for 
k = 0.01/i/Mpc, the error induced is a factor of 2. Nev- 
ertheless, there are other sources of error prohibiting the 
application of Eq. [T] in such case, such as non- negligible 
radiation and non-identical baryon and dark matter dis- 
tribution at z > 100. CMB packages such as CMBFAST 
and CAMB have already taken these complexities into 
account and shall be used to compare between high and 
low redshift observations. 



COSMOLOGICAL IMPACTS OF RELATIVISTIC 
CORRECTIONS 



In this section we discuss the impact of these relativis- 
tic corrections on a number of cosmological applications 
of LSS. We find that non- negligible impact may occur for 
searching for primordial non-Gaussianity through two- 
point galaxy clustering and N-body simulations of Gpc 
box size. 



c /local 



bga^H^g{a) 
Jb~T)6jhHl 

^" xO(l) . 



Ha 



1 



(28) 



It diverges when 6g — > 1. But this divergence is trivial, 
simply meaning that PMG of the local type does not af- 
fect the two-point clustering of 6o = 1 galaxies (Eq. 
and more detailed discussions in [2|). The induced biases 
for hg =2 and bg — 0.5 are shown in Fig. ID Existing con- 
straints on /nl'*' from two-point galaxy clustering [IJ] are 
at least one order of magnitude larger, so the bias induced 
by relativistic corrections is negligible for existing surveys 
and stage IV surveys such as LAMOST ^ and BOSS. 
However, deep full sky LSS surveys have the potential to 
constrain f^°l''^ with statistical error A/^°£^' = 0(1) [H. 
For these surveys, the relativistic correction induced bias 
should be included in the analysis. 

Terms linear in (f) in Eq. 1251 also induce bias to [i, 

are independent 



[lij . The two systematical biases in /j 
and comparable in amplitude 



local 
NL 



Large box N-body simulation 

Large simulation box size L ^ l/i~^Gpc is required 
to calculate the nonlinear matter power spectrum at 1% 
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level accuracy 17[. This is also the requirement for pri- 
mordial non-Gaussianity study. N-body simulation usu- 
ally adopts the Poisson equation V^(/) = iTrGpma^Sm to 
calculate the potential and hence misses relativistic cor- 
rections. The largest modes accessible to these simula- 
tions are k ~ 27r/L ^ 0.006/i/Mpc. N-body simulations 
typically begin at z = 100. For the largest modes, the 
density evolution to z = is then biased by more than 
20%, definitely non- negligible. 

This problem can be solved in the post-processing of 
simulation. Eq. [15] means that the simulated gravita- 
tional potential (f) is correct. So from Eq. [71 by simply 
adding a term —2a^{H/Ho)^{(j) a + 4>)/^o to the simu- 
lated 6m, we can recover the correct density evolution. 



their impact to the linear evolution of perturbations in 
ACDM. Major results are as follows. 

• We derive the exact analytical solution for the 
linear evolution of gravitational potential, veloc- 
ity and matter overdensity under the Newtonian 
gauge. It is valid at both the sub- and super- 
horizon scale, for arbitrary cosmological constant 
and curvature. 

• Relativistic corrections amplify the matter cluster- 
ing and introduce a k'^ scale dependence to the 
density fluctuation. The overdensity evolution is 
altered by 10% at fc = 0.01/i/Mpc from z = 100 to 
z = 0. 



Redshift distortion 

Redshift distortion is emerging as a promising probe 
of the large scale structure [18,]. Cosmological infor- 
mation encoded in redshift distortion is often quanti- 
fied by a single parameter /. The usual definition is 
/ = din Dm /din a. Through this definition, one can 
work out the impact of relativistic corrections straight- 
forwardly. 

Alternatively, one can define it as 



/ 



-9/aH _ V^W/aH 



(29) 



This definition does not assume a deterministic peculiar 
velocity-overdensity relation and is hence more general. 
We then have 



/ 



H'a/H +l/H^a^ da/H^ 
1 + C(fc,a) 



N 



1 + C{k,a) 



Here, /n is the / in the Newtonian limit. Relativistic 
corrections induce scale dependence in /, quantified by 
C(fc,a). 

The 7-index is often adopted to describe the structure 
growth [l^. If we straightforwardly extend the usual 
definition, we have 



j{a,k) = 



In/ 



lnil,„(a) 



(30) 



IN {a) 



ln[l + C{k,a)] 
lnf2„(a) 



In 



C{k,a) 
lnf2m(a) 



For most scales and redshifts accessible to the pre- 
cision redshift distortion measurement, C ~ (fc x 
lO'^/i^^Mpc)^^ <C 1. So relativistic corrections have neg- 
ligible impact on / and 7. 



• The gravitational potential and velocity evolution 
obtained under the Newtonian limit remaines exact 
with the presence of these relativistic corrections. 

• If these relativistic corrections are not included 
in analyzing the two-point galaxy clustering, con- 
straint on the primordial non-Gaussianity will be 
biased by Sf^l""^ - 1. 



• These relativistic corrections bias N-body simula- 
tions, which are based on the Newtonian approx- 
imation. For ^ h^^Gpc box size simulations, the 
largest mode of density fiuctuations is underesti- 
mated by - 20% from z = 100 to z = 0. We 
further show that this effect is correctable in the 
simulation post-processing. 

One question remains is whether we can extend the an- 
alytical solution to dark energy models. Under the limit 
that dark energy fiuctuation is negligible, Eq. [3 [8| & [9| 
hold. So it seems that we can repeat the derivation for 
smooth dark energy models up to Eq. 1131 which allows 
us to express the linear evolution with a double integral 
form. Unfortunately, numerically we find that the solu- 
tion obtained in this way does not reproduce the known 
(and correct) behavior at deep sub-horizon scale. The 
reason is that, dark energy fiuctuations are inevitable 
since w —1 [20|. It is true that these fluctuations 
can be small and negligible with respect to the dominant 
terms in Eq. [7l[8|&[9] However, they are comparable to 
(j) in Eq. [TSj By neglecting them at first hand in Eq. [3 [8| 
& [9l we miss non- negligible corrections from dark energy 
fiuctuations to Eq. [T51 
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